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Abstract 

The work deals with universal thermodynamics for flat FRW model, considering interacting 
n— fluid system. By introducing the state finder parameters, the feasible regions in the {r, s}— plane 
has been examined for the validity of the Generalized Second Law of Thermodynamics. A general 
prescription for the evolution of the horizon has been derived for general non-static spherically 
symmetric space time and then it is applied to FRW universe. A general restriction on the matter 
has been found for the validity of GSLT. 
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1 Introduction 

In 1974, Hawking [HE] discovered that a black hole(BH) behaves like a black body, emmitting thermal radiations 
with a temperature proportional to its surface gravity at the horizon and the related entropy is proportional 
to its horizon area [3j |4] . This entropy, temperature and the mass of the BH are related by the first law of BH 
thermodynamics (dM = TdS) [H |3] . Due to the above thermodynamics-geometry relationship i.e. 

(temperature, Entropy) ~ {Surface gravity, Area), 

There had been speculations in the literature that thermodynamic laws are in some how related to the Einstein 
equations. The real picture came in 1995 when Jacobson |6j was able to formulate Einstein equations from 
the thermodynamical laws. Subsequently, Padmanabhan [71 |H] showed the equivalence in the otherway. Most 
discussions on BH thermodynamics have been related to static BH. Subsequently, Hayward f9l fTO| have initi- 
ated the study of non-static BH (dynamic BH) by introducing the concept of trapping horizon. For spherically 
symmetric space-times he was able to write down the Einstein field equations in a form, termed as 'Unified first 
law' and first law of thermodynamics for a dynamical BH can be obtained by projecting this unified first law 
along a trapping horizon. Subsequently, a lot of works has been going on related to universal thermodynamics. 
These studies [Til HH HSl [HI [HI [ISl [171 [HI [HI [2Q] are mostly concentrated to homogeneous and isotropic 
FRW model of the universe bounded by apparent or event horizon and conditions for the validity of generalized 
second law of thermodynamics(GSLT) have been evaluated. 

In the present work, we formulate a general prescription of universal thermodynamics for FRW model start- 
ing from an interacting n— fluid system, by introducing, the state finder parameters, we derive the feasible region 
in the r, s— plane for the validity of GSLT. also formulating the evolution of the horizon, a general restriction 
on the matter has been derived and it is applied to the FRW model. 
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2 Universal Thermodynamics: A General Prescription 



We consider a homogeneous, isotropic, flat, FRW model of the universe. Suppose the universe is filled up 
with interaction n-fluids having energy densities and thermodynamic pressures are {pi,pi), i — 1,2, ...n. So the 
conservation equations are 

+ 'SH{p, + Pi) = Q„ i = 1, 2, ...1 (1) 

where Qi stands for the interaction term related to i— th fluid it may be positive or negative so that 
^"^-^ Qi = 0. Now combining the n— equations in ([1]) we have 

p + 3H{p + p) = (2) 

where p — X]r=i total energy density of the n— fluids and p — ^"^iPi, is the total thermodynamic 

pressure of the n— fluid. Thus equation ^ can be considered as the energy conservation equation of the result- 
ing single fluid having energy density p and thermodynamic pressure p. 



The Einstein field equations for the above interacting n— fluid system are given by 

3ij2 = p (3) 

and 

2H=-{p + p) (4) 

To study the thermodynamics of the universe as a thcrmodynamical system we start with the amount of 
energy crossing the horizon (Rh) during an infinitesimal interval of the time dt as 

~dE ^ A-KRlT^^m^dt ^ AtiRIH{p + p)dt (5) 

Then from the clausius relation 



-dE^dQ^ ThdSh (6) 

We have, 

+ (7) 

where dQ is the heat flux crossing the horizon during the time interval dt, Th and St are respectively the 
temperature and entropy at the horizon. 

Now the entropy {Si) of the matter of the universe inside the horizon can be related to its energy and 
pressure in the horizon by Gibb's equation 

TdSi ^ dEi +pdV (8) 

where Sj = Y^^=i the sum of the entropies of different matter components, V = |7ri?^ is the volume 
of the universe enclosed by the horizon and Ej = p.V is the total energy of the matter bounded by the horizon. 
Then using the conservation equation Q and after a little bit of calculations we obtain the variation of the 
inside matter entropy as 

^^jr [-i^RlH{p +p)+ AnRlip + p)R^] (9) 

where for equilibrium thermodynamics we chose the temperature of the inside matter same as that of the 
horizon. Thus combining ^ and ([9]) the time variation of total entropy (matter and horizon) is given by 

-^Rhip + p) (10) 



dt Th 

dt 



Thus generalized second law of thermodynamics (GSLT) (i.e. > 0) will always be true for ay one of 

the following two possibilities: 
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Figurel- Valid Regions in |r,s|-plane 
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(a) The equivalent single fluid for the interacting n— fluid system is in quintessence era and horizon is non- 
decreasing with the evolution. 

(b) The equivalent single fluid system is in phantom era and horizon is contracting with the evolution. 

It should be noted that the validity of the GSLT does not depend on any individual fluid nor on the inter- 
action among the fluids. 

Now we shall express the universal thermodynamics in terms of state finder parameters {r, s}, proposed by 
Sahni et al [211 HI] . According to them, statefinder parameters are defined as 



where g = — f§ is the deceleration parameter. 



1 (Pa , {r-l) 
and s = 



3(9 -i) 



(11) 



The parameters 'r' and 's' are dimensionless and 'r' forms the next step i the hierarchy of geometrical 
cosmological parameters after H and q [21] . It is speculated that these parameters together with coming SNAP 
observations may discriminate between different dark energy models. Further, trajectories in the {r, s} plane 
corresponding to different cosmological models demonstrate qualitatively different behaviour. 



Now using these parameters q, r, and s, the rate of change of the Hubble parameter can be written as 

3iJ2 \^ 2(r- 1) 



H^-{l + q)H' =-- 
So from the Einstein field equation (jU) we have 



1 



9s 



(12) 



p + p = m ( 1+ 



Thus variation of total entropy can be expressed in terms of {r, s} parameters as 



dSt, 



dt 



1 



2(r-l) 
9s 



Rh 



(13) 



(14) 



If we now draw the straight line 2(r — 1) -I- 9s = in the {r, s}— plane then region I (see Figl) corresponds 
to the valid region for GSLT to be satisfied provided horizon radius grows with the evolution of the universe 
otherwise region II will be the possible region for validity of GSLT. Therefore, we can classify regions in the 
{r, s}— plane where universe can be considered as a thermodynamical system. 
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Figurell-Valid Regions in |r,s)-plane 
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Further, for the comnionly used horizon in the literature, namely the apparent horizon given by 



we have 



1 ■ 3 , 



2(r-l) 
9s 



So 



dSt, 



1 



2(r-l) 



> 



dt Th [ 9s 

Thus we have the known result that GSLT is always satisfied at the apparent horizon. 

On the other hand, the event horizon (which is defined for accelerating model) is defined as 



dt 



and we have 



Re — H{Re — Ra) 



(15) 
(16) 

(17) 



(18) 



and it can not be expressed in terms of {r, s} —parameters. So in this case Re > Ra and region I or region 
II is the possible solution for the validity of GSLT. However, if we use the Ricci's length scale namely 



Rl = 2H^ +H ' 



as the boundary of the universe, then 

1 



1 



Rl - -RlH^iq + 2 - r) = IrIH^ - - - 



(19) 



(20) 



The argument behind i?L as the IR cut off in the holographic bound [531 HI] is that it corresponds to the size 
of the maximal perturbation, leading to the formation of a black hole. Thus in this case the possible part of 
the the {r, s}— plane for the validity of GSLT is classified as 
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1 \ 2(r — 1) 

Region A: - - (r - 1) ( 1 > and 1 + > 

^ ' 3s / 9s 



or 



3 , , ( 1 \ 2(r- 1) 

Region B : - - r - 1 1 < and 1 + -^^ < 

2 \ 3s/ 9s 

In figure 2 we have identified regions A and B which are the vahd region for the universal thermodynamics. 



3 Evolution of the horizon in general spherically symmetric black 
hole 

In this section, we study the growth or contraction of the horizon of a general non-static spherically symmetric 
black hole due to exchange of matter with its environment and then we extend this idea to the universal 
thermodynamics. Any spherically symmetric gravitational field can be written as 

dS^ = -/(r, t)dt^ + -y^dr^ + r^d^l (21) 

The horizon is located at r — such that f{rh,t) = = g{rh,t). To have a smooth system across the 
horizon we make the following Painleve-type coordinate transformation: 



dt^dt-.l- — - (22) 

fa 



so that metric (|2T]) becomes 



dS^ = -fdt^ + 2 J -^^^^ dtdr + dr^ + r^dnj (23) 



The incoming and out going radial null geodesies are defined by dS^ — — d9 — d(f> and so from equation 

yr^ill (24) 



dt Y g 

Now we define the evolving horizon by the condition 



= 1 (25) 



Thus if y/1 — g > 1 then out going light rays are dragged backwards towards origin. As a result we can 
define a function rh{t) such that g{r^^\t) — and is termed as evolving horizon. We can identify the mass 
function as 

2ni{r,t) = r[l- g{r,t)] (26) 

and hence the horizon 

2ni{rh{t),t)=rh{t) (27) 
Further due to spherical symmetry, the Hawking-Isreal quasi local mass function f25J namely, 

mHi{r, i) - ^ [1 - ff'^'V^rV^r] (28) 

is identified with m(r, t) defined in (|26p . So we can identify the function 'm' as the 'mass inside radius r at 
time t'. 
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Now taking total time derivative of both sides of equation p7| we have 



(29) 

^ ^ l67Trh{t) dt ^ ' 

where — 47rr^ is the surface area of the horizon and here dot and dash stands for partial derivatives 
with respect to 't' and 'r' respectively. 

Now comparing with the first law of black hole mechanics namely, 

k 

dm = —dA, (30) 
an 

the surface gravity on the horizon has the expression 

It is to be noted that the surface gravity k can also be obtained from the evolution equation of outward 
radial null vector 

rv^r = kr (32) 

where 

= |y|,yr^,0,o| (33) 

is the out going radial null vector. 
As the metric coefficient g{r,t) is positive definite for r > rh{t) so from p6| 2m{r,t) < r for r > rh{t) or 
equivalently 

1 - 2m'{rh{t), t)>0 I.e. r^g^rh, t) > (34) 

Hence the surface gravity on the horizon is positive definite. Further, using Maple(or similar software) one 
can calculate (for details see ref [?]) for any r 

So on the horizon. 



= 4w4r^(8^^^M.^^^''). (36) 

2 V 9'irh,t) 

where in deriving the last equation we have used Einstein equations G^j/ = SttT^^. So from equation (j29p 
we obtain 



8nr,T,J^FJ ' "^ ["'{^ (37) 



Further using ([ST)) the total time derivative of the mass function at the horizon gives 



dm{rh{t),t) 
dt 




= ATTrhT^,m\^ (38) 



This gives the total mass change i.e. contributions both from mass flux across the horizon and from the 
motion of the horizon itself. Thus due to inequality ([34]) we may conclude that as long as null energy condition 
(NEC) is satisfied both the horizon as well as the mass bounded by the horizon can never decrease. 

As a particular case if we consider the flat FRW model of the universe we can write the metric similar to 
(ESI as 
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dS^ = -(1 



H^R^)df - 2HRdRdt + dR^ + R^dfll 
■jj = Ra, the apparent horizon. 



(39) 



So the horizon is located at i? = Rh 
The mass term is defined as 



m{R,t) 




(40) 



and hence the mass bounded by the horizon is 



m{Rh{t),t) 



Rh{t) 



(41) 



2 



which is the usual Mishner-Sharp mass at the apparent horizon of the FRW model. Hence from ([37|) the 
evolution of the horizon is given by 



Hence from (fTl)) using (H^ we can say that validity of GSLT does not depend on the type of fluid we are 
considering- the only criteria is whether fluid satisfies the null energy condition or not. 

4 Summary and Conclusion: 

In this work we have presented a general formulation of universal thermodynamics for flat, homogeneous and 
isotropic FRW model of the universe. We have started with an interacting n— fluid system and have shown that 
the thermodynamics of the universe does not depend on any individual fluid nor on the interactions, it depends 
on the total matter density and thermodynamic pressure of the fluids i.e. for thermodynamical study we can 
take the interacting n— fluid system as a single fluid satisfying the energy conservation relation ([2|). Then we 
have introduced the state finder parameters and expressed the restrictions for the validity of GSLT in terms of 
these parameters. Also in the {r, s}— plane we have shown the feasible regions graphically. 
In section (111), we have derived the evolution of the horizon for a general non-static spherically symmetric 
gravitating system and we have found (see equation pSI) ) that horizon radius increases or decreases depending 
on whether null energy condition is satisfied or not. Then we have calculated the time derivative of the horizon 
radius (see equation |42]) for the FRW model. 

Therefore we can conclude that for a general prescription of univerasl thermodynamics in FRW model, we need 
to know the feasible region in (r, s)— parameter plane and have to examine whether the resulting single fluid 
satisfies the NEC or not. For future work, we shall try, whether similar general features for universal thermo- 
dynamics can be obtained for a general space-time model of the universe. 
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